PENTAGRAMMA MIRIFICUM AND ELLIPTIC FUNCTIONS 



(Napier, Gauss, Poncelet, Jacobi, 



Vadim Schechtman 



FAUST. Das Pentagramma macht dir Pein? 

Ei sage mir, du Sohn der Hole, 

Wenn das dich bannt, wie kamst du denn herein? 

Introduction 



In this note we give an exposition of fragments from Gauss [G] where he dis- 
covered, with the help of some work of Jacobi, a remarkable connection between 
Napier pentagons on the sphere and Poncelet pentagons on the plane. 

This gives rise to a parametrization of the variety of Napier pentagons using 
the division by 5 of elliptic functions. 

As a corollary we will find the classical five-term relation for the dilogarithm 
in a somewhat exotic disguise, cf. 6.4. 

The author gave talks on these subjects in Toulouse on January, 2011 and in 
Moscow on May, 2011. 



§1. Napier's rules 



John Napier of Merchiston (1550 - 1617) was a Scottish mathematician and 
astrologer. His main work is Mirifici Logarithmorum Canonis Descriptio (1614). 
The reader will find in [C] some historical remarks on what follows. 

1.1. For a spherical triangle with sides a,b,c and angles a,l3,'~f, 



cos c = cos a cos 6 + sin a sin b cos 7 
1 



2 

Consider a right-angled spherical triangle with sides a, b, c and angles a, 7 = 
7r/2. Let us call its parts five quantities 

T = (a, 6, 7r/2 - a, 7r/2 - c, 7r/2 - /3) 

We consider them as positioned on a circle, i.e. ordered in the cyclic order. If we 
turn them we get the parts of another right-angled triangle: 

r' = n(r) = (a', b', . . .) = (b, 7r/2 - a, 7r/2 - c, 7r/2 - /S, a) 

{Napierian rotation). Obviously — 1. 

Repeat this once more: 

g(r) := n^{T) = (a", b" , . . .) = {n/2 - a, ti/2 - c, ti/2 - f3, a, b) 

In other words, 

(a", 6", c") = 7r/2 - c, 7r/2 - a), 

this triangle is obtained by a reflection of the flrst one in the vertex ^; call it 
Gaussian reflection. 

1.2. Napier rules: 

I. sine(middle part) = product of tangents of adjacent parts: 

sin b — tan a cot a 

II. sine(middle part) = product of cosines of opposite parts: 

sin a — sin a sin c 



§2. Pentagramma Gaussianum 

2.1. Let us draw a spherical right-angled triangle which we denote P^QiP^ 
with the angles ZQ1P3P4 = ps, ZP3Q1P4 = 7r/2, /.QiP^P^i — p4 and sides P3Q1 — 
7r/2 - p5, QiPa = 7r/2 - p2, P3^'4 = Pi- 
Let us use an abbreviated notation p' — 7r/2 — p. So, the Napierian parts of 
this triangle are 

Applying the Gaussian reflexion we get the triangles 

Ti := g'-\n) = {p'i+i,p'i+4,P'i+2,P'i+5^P'i+3) 

where we treat the indices modulo 5. We denote the i-th triangle Pi+2QiPi+3, 
where ZPi+2QiPi+3 = tt/2, p,, = Pi+2^i+3- 

So we get a spherical pentagon P1P2P3P4P5; its characteristic property is 
PjPj+2 = 7r/2, cf. the picture on p. 481 of [G]. 



Set 

Gauss' notation: 



ctj = tan^ Pi 



{a,/3,-f,5,€) = (q;i,...,q;5) 

Napier rules give: 

cos Pi = sinpi+i cospi_i 
1 = cospjtanpj+2tanpi+3 

It follows: 

{■yd,5e,ea,a(3, Pj) = (sec^pi, sec^p2, sec^pa, sec^P4, sec^P5) 

Relations: 

l + a = j6, 1 + /3 = 6e, etc. (2.2. 
Out of two quantities one can build up the remaining three, for example: 

_ 1 + q; + 7 - 1 + a I + 7 

q;7 70; 

etc. and also 



7 = — , S — a/3 — 1, e 



— 1 ' ' a/? — 1 ' 

etc. One has 

3 + a + /3 + 7 + 5 + e = a/375e = ^(1 + a)(l + + 7)(1 + + e) 

2.2. Gauss' favorite example: 

{a, P, 7, (5, e) = (9, 2/3, 2, 5, 1/3), a(3-fSe = 20 

2.3. Reguleir pentagram, li a — /3 — . . . then 

q; + 1 = a^, 

whence 

1 + V5 
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a' = 11±^ = 11.0901699... 
2 

Also = sec^ i, so 



cospj = q; 



-1 



a 



1 + V5 



l + a 2 

Note that 



= — = ' " = 0.618033988... 

l + CK 2 

We have 

sin27r/5 = 2s • c 



a 



1 + V5 



where we have denoted s = sin7r/5, c = cos7r/5. On the other hand, 

sin27r/5 = sin37r/5 = s(3 - 45^) = s{4:C^ - 1), 
whence 2c = 4c^ — 1, so 

1 + \/5 a 



It follows that 
We set 

Then 



cos7r/5 

' 4 2 

2 cos Pi • cos7r/5 = 1 

c := cos27r/5 = : 

4 

cc' = 1/4, c-c' = 1/2 



A cone 

2.5. Let M be the centrum of the sphere. The vertices Pi all lie on a quadratic 
cone whith vertex M. Namely, we have 

ZPk-iMPk+i = 7r/2 

So MP3 is orthogonal to MPi and to MP5. 

We define the coordinates in such a way that M lies ar the origin, 

P3 = (1,0,0), Pi = (0,1,0) 

Then 

P5 = (0,cosp3,sinp3) 

and 

P4 = (cospi, 0, sinpi) 
Finally the ray MP2 ± P4MP5, so 

P2 = (cosp5,cosp4, -cosp3sinp5) 
The coordinates of P^ satisfy the equation 

COSP2 • (cospi • z — sinpi • x)(cosp3 • z — sinp3 • y) — xy, 

or 

2 ^ ^ 1 + q; + 7 
^ — \/axz — J^yz -— — xy :— 



— z^ -\- pxz + qyz + rxy — (C) 

Note that 

Reduction to principal axes. 



2.6. In general, given a quadratic form defined by a symmetric matrix A, 
V — {x, y, z)* I-)- v^Av, if we want to find an ortfiogonal matrix B sucfi tfiat for 
v^Bv', v'^{x',y',zr, 

v*Av = v'*B*ABv' = v'*A'v' 

such that A' = diag(G', G', G") then A' = B*AB = B-^AB and the numbers 
G,G', G" are the eigenvalues of A. 

The columns of B are eigenvectors of A. 

2.7. In our case, for the quadrics (C), the matrix is 

r/2 



r/2 q/2 
.p/2 q/2 1 



The characteristic polynomial is 

det(t ■ I -A)=t^ -t 



3 .2 + T + iij)q - r) 
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so the characteristic equation takes the form 

t(2t - If ^ uj{t - 1) {E) 

where 

This is our main equation. 

2.8. Let us investigate the real roots of (E). We suppose that a; > 0. The 
straight line 

£: u^uj{t-l) (2.8.1) 

allways intersects the cubic parabola 

P : u^t{2t-iy (2.8.2) 

at one negative point t = G < 0. 

On the other hand, i intersects P at two points t = G', G" > 1 iff a; is greater 
than some critical value ujq] ii ui — ujq then I touches P &X, t — G' — G" , ii ui < ujq 
then there are no points of intersection of £ with P other than t — G. 

The critical value is ouo — ckq where 

1 + ^5 

«o = ^— , 

cf. 2.3. In that case 

G — —ceo, Gf' — G" — q;q/2 = — cqG, 
Co = cos7r/5 



Thus, if a; > q;q then {E) has one negative root, G, and two positive roots, G', 
G" which coincide \l u — a^. 

One has: 

GG'G" = -u/A 
(G'-l)(G"-l)(G'"-l) = -l/4 
(2(7 - 1)(2G" - 1){2G" -l)^-uj 

2.9. Example. For a^-fSe = 20, G = -2.197, G' = 1.069, G" = 2.128. 



§3. Gauss' coordinates 



3.1. Consider an eUipse 



^4 + S = i. 



a > b. Inscribe E into the circle C with the centrum O and radius a. For a 
point P = {x,y) G ii^, let P' = {x,y') G C. The eccentric anomaly (anomalie 
excentrique) of P is the angle = ZXOP' where X = (a, 0). Then 

X = a cos (p, y — b sin 

3.2. Return to the Napier pentagon Pi . . . P5. Let us draw a plane tangent 
to the sphere at the point of intersection with the axe of the cone and take the 
central projection of our pentagon on this plane. Let Ri be the projection of Pi to 
this plane. The points Ri, i = 1, . . . , 5, he on an ellipse with semi-axes ^/—G/G' 
and ^J-GjG". 

We put the coordinate axes x, y of the plane along the axes of the ellipse. Let 
Ri have the coordinates (xj, yi). Let 0^ be the eccentric anomaly of Ri. We have 

Xi = g' cos 0i, yi = g" sin 0j (3.2.1) 

where 



Let M be the centrum of the sphere, with coordinates {x, y, z) = (0, 0, 1). 
Attention: these coordinates differ from 2.5. 
Set = ZPiMPi+i = ZRiMRi+i] = (xi,yi, 1). 
Then 

(ri, ri+i) XiXi+i + yiyi+i + 1 

cos a,- = — 



ml \/(^' + l/' + l)(a^^i + 2/m + l) 



and 

ai :— tan^ ■^j = sec^ ipi — ^ 



For the future use, note the quantities 



(3.2.2) 



Ui + l |rj|^|rj+i|^ 

3.3. We have 
whence 

Xi-iXi+i + yi-iUi+i + 1 = (3.3. l)i 

Solving (3.3.1)j_i and (3.3.1)j+i for Xi,yi, we get 

yi+2 — yi-2 — /„ n o\ 

Xi = , yi = (3.3.2) 

Xi+2yi-2 - yi+2Xi-2 Xi+2yi-2 - yi+2Xi-2 

3.4. Next, 

XiXi^^ _mh+^ I ^ , . 

2G'-1 2G" -I 2G-1 ^ ' 

This will be proven in 3.6. It follows: 

2G' - 1 i/i+i - yi_i 



2G - 1 Xi-iyi+i - Xi+iyi^i ' 



2G" - 1 Xi+i - Xi^i 

yi = -7^ — 3-4-2) 

2G - 1 Xi-iyi+i - Xi+iyi-i 

3.5. Theorem (Gauss, April 20, 1843). (a) 

sin((0,_2 + 0.+2)/2) _ G _ cos((4-2 + 0.+2)/2) _ G 

cos((0,_, - 0,+2)/2) " G" cos((0,_, - 0,+2)/2) ~ G' '^''^^ ^'-'-'^ 

(b) 



sin((0,_i + 0,+i)/2) _ / G{G-l) . _ G{2G - 1) . 
cos((0,_i - 0,+O/2) " V G"{G" - 1) '^^'^^ - G"{2G" - 1) '^"'^^ 



cos((0.-. + 0..O/2) ^ ^^(^-1) ^^,^^ ^ G(2G-1) 

cos(((/),_i - 0,+i)/2) YG'(G'-l) G'{2G'-l) ^ ^ 

3.6. Proof, cf. [F]. (a) follows from (3.3.2) and (3.2.1), taking into account 
the formulas 

sin a — sin 6 = 2cos((a + h)/2) sin((a — b)/2 
cos a — cos 6 = — 2sin((a + h)/2) sin((a — h)/2 (3.6.1) 

Proo/o/ (3.4.1). 



It follows from (3.5.1) (replacing i + 2 by i) that 

G'^ cos\{4>i + 0i+i)/2) + G"^ sin2((0, + 0,+i)/2) = G^ cos\{4>i - 0i+i)/2). 
Using 

2 1 + cos 2a .2 1 ~ cos 2a 

cos a — , sm a — 

2 2 

we get 

cos(0, - 0,+i) + (C'^ - G"^) cos(0, + 0,+i) = G"^ + - 
Prom (3.2.1) we get 

_ ^ G"2)G'a;,x.+i + (G^ + G'^ - G"2)G"|/,y,+i+ 

+(-G2 + G'" + G"')G = (3.6.2) 
The numbers G, G', G" are the roots of 

t{2t - If = u{t - I) (3.6.3) 

It follows that 

G^^G'^ + G"^^^^ 

Whence 



2 2G-1 

Similarly we compute the coefficients at XjXj+i and at yiUi+i of (3.6.2) and arrive 
at (3.4.1). 

Proof of (b). The second equalities follow from (3.6.3). The first ones follow 
from (3.4.2) in the same manner as one has deduced (a) from (3.3.2). □ 



§4. Poncelet's problem and division of elliptic functions (Jacobi) 



In this Section we describe the work of Jacobi [J2] . For a modern treatment 
of it see [GH] and references therein. 

4.1. One considers a circle with center G of radius R, inside it a smaller circle 

with center c of radius r, a will be the distance between their centra. The line 
cG intersects the bigger circle at a point P, so \CP\ = R, \cP\ = R + a. 

One takes a point Ai on the bigger circle, draws the tangent to the smaller one 
till the intersection with the bigger one at the point ^425 and continues similarly. 
Denote 

2(f)i = ZAiCP 



Let us suppose for simplicity that C is inside the smaller circle. Let B be the point 
where A1A2 touches the smaller circle, B' e A-1A2 the base of the perpendicular 
from C, and D e cB the base ofthe perpendicular from C. Then 

= CS' = i? cos(02-0i) 

and 

cD — acos(02 + 0i), 

so 

r = cD + DB = Rcos{4>2 — 4>i) + acos(02 + 4>i) 

It follows that 

it:cos(0j+i - 0i) + acos(0j+i + 0j) = r 

or 

{R + a) cos cos (t)i-\- {R — a) sin sin (f)i — r (4.1.1)j 
Substracting (4.1)^ — (4.1)i_i and using 

cos X — cos y 



sm y — sm X 

we get 



tan 



/ x + y \ 
V 2 J' 



tan((0i+i + 0i-i)/2)) = ^— tan^i (4.L2) 

4.2. Jacobi elliptic functions. Cf. [Jl]. Fix a number < A; < 1 and define 
the period 

Jo ^/l-k^sm^x Jo ^{1 - y^){l - k^y^) 
The function "amplitude" 

am : [0,K] — > [0,7r/2] 

is defined by (f> — am(M), 

dx r^'f' dy 



u 



lo Vl-k^sin'x Jo V(l-2/')(l-^V) 
Thus 

am(X) = 7r/2 

We extend am to a function R — > R by a,m{—u) — —a,m(u), a,m{u + 2K) 
am(u) + TT. 

We set 

(iamw /- . 2 

A am 14 := — ; — = v 1 — A;"^ sm am u 
du 

Modern notations: 



sn = sin am 14, cnii = cosamw, dnw = Aamw 
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Thus sn is an odd function and cn and dn are even. 
4.3. Addition formulas. 



snucnvdnv + cnusnvdnu 

sn{u + v) = — 

1 — sn^ u sn^ v 

cn M cn f — sn -u sn f dn M dn 1) 



cn{u + v) 
dn{u + v)- 



1 — sn^ u sn^ v 
dn M dn f — /c^ sn m sn f cn -u cn v 



1 — k'^ sn^ u sn^ v ' 
cf. [Jl], §18. 

Main formula: 

cn(M — v) = ciaucnv + snusnv dn{u — v) 

Next, 

2 cn -u sn f dn -u 
sii{u + v) - sn{u -v)^ — — — 

1 — k^ sn"' u sn^ v 

2 sn sn i> dn 14 dn 
cn(M + v)- cn{u - v) ^ — — — , 

1 — k^ sn^ u sn^ v 

whence 

// / N ^ \\ cnx — CUT/ 

tanffamfx) + am(y))/2) = = 

sny — snx 

dn{(x — y)/2) tanam((x + y)/2) 
4.4. If (pn = am(0 + nt) then 

tan((0o + 02)/2) = A am t tan 01 

So if 

R-a 

A am t — — 

R + a 

then (4.1.2) is satisfied. 

We also have 

cos (pi cos (pi+i + sin (pi sin • y^l — k^ sin^ a — cos 
where a — amt. 

So we can find k and a from the equations 

^/l — k^ sin^ a — — 

R + a 

and 
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wherefrom 

k-^ ^ (^-af-r^ 

(i? + a)2-r2 (i? + a)2-r2 ^ " " ^ 

It follows: 

4.5. Theorem. // the process closes up after n steps and m full turns then, 
defining k through (4.4.1), we have 

dd) m 



•\/r^-^^"sii?0 ^ Jo \Jl-k'^ sin^ 



§5. Back to pentagramma 



5.1. Let us return to (3.3.1). Denoting = ^j, 0' = 0j_|_2 and using (3.2.1) we 
get 

g'^ cos cos 0' + g""^ sin sin cf)' — 1 

or 

COS0COS0' + — sin0sin0' = — — (5.1.1) 

Gr Gr 

Now compare this with the Main formula (4.3.2). We see that (5.1.1) will be 
satisfied ii (f) — eim{u), (f)' — am{u + w), a — am{w), 

Q' 

dnw = V 1 - /c^sin^Q! = — (5.1.2) 



and 



It follows that 



and 



G' 

cos a — ~7T- 



G' 

cnw ^ -— (5.1.3) 
G 



_ rill -2 

We can go backwards. 

5.2. Theorem. For < k < 1, let K the corresponding complete elliptic 
integral (4.2.1). Define vectors in R^; 

r(ku)= /^£^(!i±iW^ Vdn(2X/5) sn(^ + 4jX/5) \ 
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M G M, j G Z. Then rj{k,u) = Yj+^{k:,u). Set 

_ \Yj{k,u) X rj+i(/c,M)p 

T/ien u;e have 

1 + Q;j(A;, li) = q;j_2(A;, u)aj+2(A;, li) 



(5.2.1) 



In the degenerate case /c = we have 

for all j, u. The ends of vectors rj(0, li) form a regular plane pentagon, cf. 2.3. 

5.3. The idea of uniformising relations in the spherical (resp. hyperbolic) 
geometry by elliptic functions (cf. [S] for a nice review) plays an important role 
in Statistical Physics. Onsager uses it in his 1944 paper [O] on the solution of 
the Ising model. Baxter remarks that the same idea underlies the solution of his 
star-triangle relations, cf. [B], 7.13. 

Baxter cites a book [Gr] in this connection; apparently this is the very same 
" Greenhill's very odd and individual Elliptic functions^'' which, according to Lit- 
tlewood cited by Hardy, was the Ramanujan's textbook, cf. [H], Ch. XII. 



§6. Dilogcirithm five-term relation 



6.1. Euler dilogarithm: 

r log(l-t) 



oo 



n=l 



0<x<l. 

Rogers dilogarithm: 



< X < 1. 

6.2. Theorem (W.Spence, 1809). 



L(x) + L(l - x) = y (6.2.1) 

L{x) + L{y) - L{xy) - L{x{l - y)l{l -xy))- L{y{l - x))/(l -xy))= 0, (6.2.2) 
< x,y < 1. 
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6.3. One can reformulate (6.2.2) as follows (cf. [GT]): 

L{x) + L(l - xy) + L{y) + L((l - y)/(l - xy)) + L((l - x)/{l - xy)) = 
Define 

. . . , 65) = (x, 1 - xy, y, (1 - 7/)/(l - xy), (1 - x)/(l - xy)) 
and for an arbitrary n e Z define 6„ by periodicity bn = bn+5- 
Then 

= 1 - n e Z. 

Define a new sequence a„ by 



1 - 6n ' 1 + On 

Then 

Combining this with 5.2 and (3.2.3) we get 
6.4. Corollciry. Let rj{k,u) be as in 5.2. Set 

\rj{k,u) X r^+i(A;,M)p 



I3j{k,u) 
Then 



\rj{k,u)\'^\rj+i{k,u)[^ 



j^Hmu))^^ (6.4.1) 



In the regular case A; = (6.4.1) takes the form 



Lia-') = - (6.4.2) 



(Landen), cf. [K]. Here as usually a = (1 + \/5)/2. 
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